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Where is the 3 subgroup of GRBs? 
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Abstract. It is shown that in the duration-hardness plane the GRBs 
of the third intermediate subgroup are weh defined. Their durations 
are intermediate (i.e. roughly between 2 and 10 seconds), but their 
hardnesses are the smallest. They are even softer than the long bursts. 
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1. Introduction 

It is widely accepted that the short and long gamma-ray bursts 
(GRBs) are really different phenomena (see, for example Norris et al. 
(2001) and the references therein). In 1998 Horvath made a trimodal 
fit of the BATSE Catalog bursts' duration distribution and found a 
third subclass of GRBs' (Horvath 1998). Later papers (Mukherjee et 
al. 1998, Hakkila et al. 2000, Rajaniemi & Mahonen 2002, Horvath 
2002) used more parameters (e.g. peak fluxes, fluences, hardness 
ratios), and in this multidimensional space studies all of them con- 
firmed that the third bursts population was statistically necessary. 

Bagoly et al. (1998) showed that in this high dimensional pa- 
rameter spaces only two main parameters were necessary to char- 
acterize all the BATSE Catalog bursts' properties. Hence, a two 
dimensional space looks like a good characterisation of the GRB 
subgroups. Therefore, here we use Tgo and the hardness H32 ratio 
in our newest analysis. Figure 1. shows the observed BATSE bursts' 
distribution on the log(T9o) — log(i?32) plane. 
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2. The fits 

This distribution can be fit- 
ted at the first step by the sum of 
two two-dimensional normal dis- 
tributions in the logTgo - logii"32 
plane. If one were fitting si- 
multaneously the values of logTgo 
and logi?32 by one single two- 
dimensional (bivariate) normal 
distribution, then the distribu- 
tion would have five independent 
parameters (two means, two dis- 
persions, and the correlation co- 
efficient). The standard form of such a bivariate distribution is given 

by 

f{x,y)dxdy = 
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Fig. 1. The values of logTgo and 
logi732 of the 1929 BATSE GRBs. 
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where x = logTgo, y = logi?32, C = 2r{x - ax){y 
means, ax, o-y are the dispersions, and r is the correlation coefficient. 
N is the number of GRBs, and f(x,y)dxdy the theoretically expected 
numbers of GRBs at the inf. surface at the [x,y] plane given by 
intervals [x,{x + dx)] and [y, (y + dy)]. In other words, {f (x , y) / N)dxdy 
defines the probability of finding a GRB at the given inf. surface. 

This fit was done by a stan- 
dard search for 11 parameters with 
N = 1929 measured points. Each 
GRB defines a point in the x,y 
plane with coordinates Xi,yi {i = 
1,2,...,N). The theoretical curve 
f2ix,y,axk,ayk,(Jxk,cryk,rk,W) (k = 
1, 2) is a sum of two normal distribu- 
tions as given in Eq.l. The normal- 
ization constant of the first [second] 
term is NW [N{l-W)], where W is 
the weight of the first normal distribution. For the first (second) 
term the parameters are axi,ayi,axi,cryi,ri (a^2 , ay2 , o-x2 ,o-y2,r2)- 



Table 1. Best fit with two 
bivariate normal distributions 
for logTgo a-i^d logi/32. 
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We obtain the best fit to the 11 parameters through a maxi- 
mum hkehhood (ML) estimation. We search for the maximum of 
the formula 
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using a simplex numerical procedure; the index "2" in L2 shows that 
we have a sum of two log-normal distributions of type given by Eq.l. 

The results of this fit are shown in Table 1. One can calculate 
a density distribution of the observed data on the log(T9o)-log(H32) 
plane. The values of Table 1 define the theoretical distribution of 
GRBs, if there is any theory, which suggests lognormal H32 and T90 
distributions, in the Tgo - i?32 plane under the assumption that there 
are only two subgroups. 

Hence, one can make the dif- 
ference between the actual distri- 
bution and the theoretical distri- 
bution. If there were only two 
components, which distributed 
lognormally, one would find only 
some noise being proportional to 
the square root of the actual 
value of the burst density func- 
tion. Surprisingly the biggest de- 
viation is not there where the ob- 
served (or the theoretical) burst 
density function is the biggest, 
not even the dense core of the two 
subgroup population, rather then 
between the two groups in T90 but 
in H32 softer then either groups. 

This deviation is shown by 
Figure 2. The widest contour has the half intensity then the highest 
deviation. In order to be sure that this difi'erence is not a chance and 
is not given by the random noise, we proceeded as follows. We take 
this part of the logTgo - logi?32 plane. This part contains 145 bursts. 
The integral over this area of the theoretical distribution gives 63. 

Depends on which number is the Poisson parameter one can 
calculate the probability 10~^'''-10"^^. Of course there are so many 
different area and shapes, which we can probe. However if we take 
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Fig. 1. Departure of the numbers 
of GRBs from the theoretical values 
given by the best fit collected at Ta- 
ble 1. 
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10 different shapes 100 different sizes and try 1000 different positions 
in tlie plane tlie probability is still very low. Therefore one can say 
the deviation which the Figure 2. shows is very unlikely can cause 
by chance. 

Unfortunately, the possible third group members are mixed with 
the long ones (partially also with the short ones, too). 

3. Conclusions 

We have argued that both the duration and also the hardness 
should be distributed log-normally - of course, in any subclass sep- 
arately. We also provided in the logTgo - logiJ32 plane a fit with 
the sum of two bivariate normal distributions. The idea for this 
fitting was given by the observational fact that the short and long 
subclasses are different both at hardnesses and durations. Finally 
we obtained the difference between the actual distribution of GRBs 
and the theoretical one. The difference, not being a Poissonian noise, 
shows another subgroup, the third subgroup. 

As the result we obtained the intermediate subclass having ap- 
proximately the same number of bursts than it was previously pre- 
dicted ( Horvath 1998, Mukherjee et al. 1998, Hakkila et al. 2000, 
Rajanicmi & Mahonen 2002 and Horvath 2002). We also confirm 
that their hardnesses are low. 
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